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ABSTRACT 



on . 

Aims. The statistical equilibrium of neutral and ionized silicon in the atmospheres of metal-poor stars is discussed. Non-local thermodynamic 
O ' equilibrium effects are investigated and the silicon abundances in metal-poor stars determined. 

Methods. We have used high resolution, high signal to noise ratio spectra from the UVES spectragraph at the ESO VLT telescope. Line 
formation calculations of Sii and Sin in the atmospheres of metal-poor stars are presented for atomic models of silicon including 174 terms 
and 1132 line transitions. Recent improved calculations of Sii and Sin photoionization cross-sections are taken into account, and the influence 
of the free-free quasi-molecular absorption in the Lya- wing is investigated by comparing theoretical and observed fluxes of metal-poor stars. 
All abundance results are derived from LTE and NLTE statistical equilibrium calculations and spectrum synthesis methods. 
Results. It is found that the extreme ultraviolet radiation is very important for metal-poor stars, especially for the high temperature, very metal- 
poor stars. The radiative bound-free cross-sections also play a very important role for these stars. 

Conclusions. NLTE effects for Si are found to be important for metal-poor stars, in particular for warm metal-poor stars. It is found that these 
effects depend on the temperature. For warm metal-poor stars, the NLTE abundance correction reaches ~ 0.2 dex relative to standard LTE 
calculations. Our results indicate that Si is overabundant for metal-poor stars. 
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. ^ ' 1. Introduction (LTE). They show pronounced star-to-star scatter at all metal- 
s' . licities, and particularly for extremely metal-poor stars, where 

. Silicon is not only an important reference element for compar- u jo-uj • .. ju i <, r a 

H the observed Si abundance is represented by only two lines and 

d ' ision of various types of cosmic matter with the Sun, but also ff , ■ , . ... e PIJ , u ,. r™ - • 

. . . JF affected by the contamination of CH and Ha lines. This is seen 

one of the main electron contributors and opacity sources in the ■ „• m„ . . , i onri£ i , T • . , ^ mo i e A tU t 

r ■> in FigUJ Preston et al. (2006) and Lai et al. (2008) found that 

near UV in the atmospheres of cool stars. With few exceptions, , ■ , c . u , , , , ^ j. i 

F F ' the derived Si abundance depends only on / e ff and not on log 

silicon also represents the only or-element measurement^ a ft and Preston gt a , (2006) suggested that the si abundances 

given damped Lyo* system (DLA, Prochaska & Wolfe EO®. for ^ ^ ^ > 5g(X) R probably do QQt represent the ttue 

Although silicon is a refractory element, its depletion is mild in yalues Cohen ^ a , noted ^ ^ observational picture 

lightly depleted regions of the ISM, and therefore the observed t c - ■ , , ■ r o-rc i c ■ * 

6 3 F & for Si shows a puzzle - the abundance ratio [Si/Fe] of giants 

column densities for the DLA system should nearly reflect the , . n . , , . ■ , t , e , , 

3 3 is about 0.4 dex larger than that obtained for dwarfs (also see 

total Si column density. Ryan et al. 119961 ). and they discussed in detail the possibility of 

The current paradigm is that silicon is made during oxygen tb j s difference - 

and neon burning in massive stars, and it would therefore be 

ejected by Type II supernove (Woosley & Weaver [l995l Umeda - There is a problem for the transition probabilities of the two 
et al. 120551 Ohkubo et al. [2pQ6l . Type la supernove also pro- very strong Si i lines. Cohen et al. (|2004|) studied high tem- 
duce some sihcon (Tsujimoto et al. [T995l Iwamoto et al. [l999l perature extremely metal-poor dwarfs, but they used only 
Previous abundance determinations of Si have been carried the 3905 ° K Si 1 line - This line overlaps a CH f eature ; how- 
out under the assumption of local thermodynamic equilibrium ever > 11 is ver Y weak in such stal " s - Ca y rel e J al - 12QQ4]) relied 
instead on a single line of Sii at 4102.9 A, which is close 
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Fig. 1. Silicon abundances of metal-poor stars determined from 
LTE analyses of Fulbright (2000, Filled circles • for dwarfs 
and open circles o for giants), Mc William et al. (119951 pluses 
+), Honda et al. ( 120041 crosses x), Aoki et al. ( 120051 asterisks 
*), Ryan et al. ( 119961 filled diamonds ♦ for dwarfs and open 
diamonds o for giants), Gratton et al. (2003, filled squares ■), 
Cayrel et al. ( 120041 open triangles A), and Cohen et al. (120041 
filled triangles a). 



value from the laboratory work of Garz (119731) . However, 
NIST gives -2.92 dex for the 4102 A line, which is 0.22 dex 
higher than the laboratory value. Thus, the Si abundance 
derived by Cayrel et al. (2004) would be 0.22 dex lower if 
the NIST gf value were adopted. 
- Only one non-local thermodynamic equilibrium (NLTE) 
analysis has been published to date, in which Wedemeyer 
(2001) reported the influence of deviations from LTE upon 
the Si abundances in the Sun and Vega, with no reference 
to line formation in metal-poor stars; he demonstrated that 
NLTE for Si in the Sun is negligibly small. However, he did 
not consider line formation for the two strong lines. 



Cohen et al. (2004) suggested that this difference probably 
arises from contamination of the Si i line at 3905 A by blend- 
ing CH features in the giants; spectral syntheses were not car- 
ried out for the giants by Cayrel et al. (2004). Preston et al. 
(2006 ) examined this suggestion in their red horizontal -branch 
(RHB) stars and found that the CH band is very weak. Testing 
synthetic spectra of cooler RHB stars computed using derived 
carbon abundances revealed that the CH contribution can only 
be a few percent of the total absorption at the Si i feature. They 
also found that the derived Si abundances for cooler stars are in 
reasonable accord with values determined for stars of similar 
metallicity using Si i transitions in the red spectral region. 

The present work is based on a sample of metal-poor stars 
and aims at exploring their [Si/Fe] abundance ratios applying 
a full spectrum synthesis based on level populations calculated 
from the statistical equilibrium equations. In Sect. 2 we present 
the observational technique and the atmospheric models and 
stellar parameters are discussed in Sect. 3. NLTE line formation 
is discussed in Sect. 4. The results and comparison with other 



works are shown in Sec. 5. The discussion is presented in Sect. 
6, and the conclusions are given in Sect. 7. 



2. Observations 

Our approach to obtain a representative silicon abundance in- 
vestigation involves analyzing a sample of metal-poor stars 
with high resolution and high signal to noise ratio spectra. 
The spectra of 13 metal-poor dwarf stars were observed in 
March 2001 with the UVES echelle spectrograph mounted at 
the ESO VLT2. UVES observations cover a spectral range be- 
tween 3300 and 6650 A with gaps of -100 A around 4570 A 
due to the beam splitter and near 5580 A at the edge of the but- 
ted CCD. The spectral resolution is around R = 60000. These 
observations were originally intended to show a high signal 
in the blue (see Mashonkina et al. 2003). Consequently, the 
green/red spectra have a S/N near 300 in most of the single 
exposures. We also use high-quality observed spectra from the 
ESO UVESPOP survey (Bagnulo et al. [20031 for HD 122563. 
For G64-12 we use the spectrum from the High Dispersion 
Spectrograph at the Nasmyth focus of the Subaru 8.2 m tele- 
scope (Noguchi et al. 120021) . 

Data extraction followed the standard automatic IDL pro- 
gram environment designed for the FOCES spectrograph 
(Pfeiffer et al. 1998), but with slight modifications also ap- 
plicable to the UVES data. All echelle images including flat 
field and ThAr were corrected for bias and scattered light back- 
ground. Objects and ThAr exposures were extracted and cor- 
rected for flat field response. Bad pixels were detected and as 
far as possible removed by comparison of the 3 single expo- 
sures (see Gehren et al. 2004 for detail). 



3. Method of calculation 

3.1. Model atmospheres 

Our analyses are all based on the same type of atmospheric 
model, irrespective of temperature, gravity or metal abundance. 
We use line-blanketed LTE model atmospheres, generated as 
discussed by Fuhrmann et al. ( 11997I ). The main characteristics 
are: the iron opacity was calculated with the improved mete- 
orite value log e Fe = 7.51 (Anders & Grevesse 1989); opacities 
for metal-poor stars with [Fe/H] < -0.6 were calculated us- 
ing ff-element (O, Mg, Si and Ca) abundances enhanced by 0.4 
dex, and the mixing-length parameter l/H p was adopted to be 
0.5, in order to determine consistent temperatures for H ff and 
the higher Balmer lines (see Fuhrmann et al. 1993 ). 



3.2. Stellar parameters 

For most of our program stars, we adopted the stellar param- 
eters determined by Gehren et al. (2004, 2006), where the ef- 
fective temperatures are derived from the wings of the Balmer 
lines; surface gravities are based on the HIPPARCOS paral- 
laxes. Iron abundances are obtained from Fe n lines, and the 
microturbulence velocities are estimated by requiring that the 
iron abundance derived from Fe n lines should not depend on 
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Table 1. Atomic data of silicon lines*. 
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-29.869 
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-5/'D 2 
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4? 2 S 
As S 1/2 


-4p 2 P° /2 
"4p 2 P 1/2 


0.26 
-0.06 


-30.200 
-30.200 



"logg/ values have been determined from solar spectrum fits, and 
damping constants log Ct, for Si i lines are computed according to the 
Anstee & O'Mara (1991, 1995) interpolation tables. 



equivalent width. For HD 84937 and HD 122563 the parame- 
ters were taken from Mashonkina et al. ( 120081 1. where the stel- 
lar effective temperatures were determined from the hydrogen 
H a and Hy3 line wing fitting based on NLTE line formation; 
surface gravities are based on the HIPPARCOS parallaxes, and 
iron abundances are derived from Fe n lines and the microtur- 
bulence velocity was derived from the strongest Fe n, Ca i, and 
Mgi lines. For consistency with the whole series of our ear- 
lier NLTE studies (from Baumiiller & Gehren [T997l to the most 
recent paper of Shi et al. 120081) . we adopt effective tempera- 
tures based on the hydrogen resonance broadening calculated 
with Ali & Griem (1966) theory. We note that two papers of 
Barklem et al. 120001 and Allard et al. dSOOgi l obtain larger val- 
ues of the hydrogen self-broadening cross-sections. An impact 
of using the self-broadening formalism of Barklem et al. (2000 ) 
on effective temperature determinations for metal-poor stars is 
discussed by Mashonkina et al. (2008). The uncertainties for 
the temperature, surface gravity, metal abundance and micro- 
turbulence velocities are generally assumed to be +50 K, 0.05 
dex, 0.05 dex and 0.1 km s _1 respectively. 

3.3. Atomic line data 

Table [T] lists the relevant line data with their final solar fit val- 
ues (Shi et al. 2008 hereafter paper I). Collisional broadening 
through van der Waals interaction with hydrogen atoms is im- 
portant for strong Si lines. As already pointed out by Gehren 
et al. (120011 [2004), the resulting values of the van der Waals 
damping constants are mostly near those calculated according 
to Anstee & O'Mara's (fTWTl [19951 ) tables. In our analysis, 
the absolute value of the oscillator strengths is unimportant be- 
cause the abundances are evaluated in a fully differential way 
with respect to the sun. 



4. NLTE calculations 

4.1. Atomic model 

The silicon model atom includes the most important levels of 
Si i and Sin and comprises 132 terms of Sii, 41 terms of Sin, 
plus the Sim ground state. The atomic properties are docu- 
mented in Shi et al. (paper I). Oscillator strengths and photoion- 
ization cross-sections were calculated by Nahar & Pradhan 
(119931 ). The full analysis of the solar spectrum (Kurucz et al. 
1984 ) allows a reasonable choice of the hydrogen collision en- 
hancement factor resulting in Sh — 0. 1 for Si (see paper I). 

All calculations are carried out with a revised version of the 
DETAIL program (Butler & Giddings 1985) using accelerated 
lambda iteration (see Gehren et al. 200T1 120041 for details). 

4.2. NLTE effects 

The abundance analyses of silicon clearly show the NLTE ef- 
fect. There is a tendency that the NLTE effect is large for warm 
metal-poor stars, as would be expected. Our results confirm 
that LTE abundances can be significantly different from their 
NLTE counterparts, with differences reaching more than 0.20 
dex in extreme cases. It is evident that the NLTE effects are sys- 
tematically stronger for hotter models, which is in agreement 
with the statistical equilibrium of sodium, magnesium and alu- 
minium (Shi et al. 120041 Zhao & Gehren [20001 Baumiiller et 
al. |1997l ). As expected, the strongest departures from LTE are 
found for models with high temperature and low metallicity. 
For our program stars, the differences of abundance between 
LTE and NLTE analyses for the strong 3905 and 4102 A lines 
are plotted in Fig. [2] as a function of metal abundance, temper- 
ature and surface gravity, respectively. 

It should be noted that the NLTE effects differ from line 
to line, reflecting their individual properties. For example, the 
strong 3905 and 4102 A lines show large NLTE effects com- 
pared to the other weak lines. From Table [2l we can see the 
importance of the NLTE effects on different lines in the silicon 
abundance determination: 

- the weak lines show the smallest NLTE abundance effects 
(< 0.02 dex) for the stars studied here, so they are the 
best abundance indicators for LTE analyses of moderately 
metal-poor stars. 

- in the 3905 A line, which is very important for determining 
the silicon abundances of extremely metal-poor stars, the 
NLTE correction is relatively small (~ 0.1 dex) for moder- 
ately metal-poor dwarfs; it increases for extremely metal- 
poor warm stars (> 0.15 dex). Therefore, this line should 
no longer be analyzed with LTE. 

- the 4102 A line is also important in determining the silicon 
abundances of metal-poor stars; it is necessary to take ac- 
count of the NLTE effects, except for very cool and metal- 
rich stars. For warm metal-poor stars, the NLTE correction 
can reach more than 0.25 dex. 

- for the two Sin lines, NLTE abundance corrections are 
negative. A similar behaviour is found for Can lines 
(Mashonkina et al. 120071 ). 
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Fig. 2. Difference of [Si/Fe] abundance ratios calculated under 
NLTE and LTE assumptions for the two strong Si i 3905 and 
4102 A lines as a function of metal abundance (a), temperature 
(b), and surface gravity (c). 



Our results show that departures from LTE of the silicon 
level populations appear to be larger for warm metal-poor stars. 
Also, there is a clear trend that the NLTE effects increase with 
increasing temperature, which can partly explain the observed 



Fig. 3. The UV flux emerging from MAFAGS-ODF model for 
HD 84937 without (dashed line) and with (solid line) the Lya 
continuum absorption. 



anomalous T e ff dependence of [Si/Fe] (Preston et al. 2006} Lai 
et al.|2008). 

4.3. The influence of the quasi-molecular absorption 
in Lya 

Recent calculations using the R-matrix code provide detailed 
photoionisation cross-sections with autoionization resonances 
for the most important energy levels of Si i and Si n. For Si i, the 
cross-sections of the lowest three terms, 3p 3 P, 3p l D and 3p l S 
are an order of magnitude greater than those of the next terms. 
It tends to efficiently decouple the metastable terms from the 
excited ones. The thresholds of the three terms are 1500, 1700 
and 2000 A (Paper I). In this wavelength region, the most im- 
portant continuous absorption is the free-free quasi-molecular 
absorption and satellites in Lya due to collosions with H and 
H + (Allard et al. 1998), especially for warm metal-poor stars. 
In FigfJ] we compare the synthetic ultraviolet (UV) fluxes for 
HD 84937 without and with this process considered. When this 
effect is included, the flux near 1700 A decreases by more than 
an order of magnitude, and it allows us to reproduce the most 
characteristic UV features of warm metal-poor stars, as illus- 
trated by a comparison with IUE data (Cacciari |19851 l for the 
hot metal-poor star HD 84937 (Fig]4ji. 

As already discussed in paper I, there is an unusually large 
energy gap between the ground state, 3p 3 P, the two metastable 
terms 3p 1 D and 3p 1 S , and the first excited levels of neutral 
silicon. These gaps are ~ 5, 4, and 3 eV, respectively, and they 
shift all lines emerging from those levels into the UV, most of 
them below 2000 A, Therefore, the interaction of the three most 
populated Si i levels with the Si n ion is completely based on 
bound-free processes, where photoionisation and ionisation by 
electron collisions compete in strength. Due to the large pho- 
toinisation cross-sections for the three lowest levels, the depar- 
ture coefficients of Sii are very sensitive to the UV radiation 
field. When the Lya continuum absorption is not included, the 
extreme underpopulation of Si i levels in the atmospheres of 
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Fig. 4. The comparision of the UV flux from MAFAGS-ODF 
model (solid line) with the IUE observations (dashed line) for 
HD 84937. 

warm metal-poor stars is clearly seen. The depopulation with 
respect to the LTE case already starts in very deep photospheric 
layers, with a very large net pumping to the Si n ground state 
due to the strong UV radiation. For Sh = 0.1 the departure 
coefficients b, drop near to at log t = 0.5 (Fig|5] top). As ex- 
pected, the 3p l D and 3p l S terms follow 3/? 3 P tightly. This 
would not affect neutral Si line formation very much, except 
that the strong 4102 A line turns into emission for warm metal- 
poor stars, such as HD 84937, as shown in Fig|6] Also, the 
NLTE correction would be very large (~ 1.0 dex) for the 3905 
A line. 

Models with Lyor continuum absorption included produce 
much smaller departures (Fig|5] bottom), as radiative processes 
(photoionisation) are compensated by hydrogen and electron 
collisions. Since all the excited terms are only loosely coupled 
to each other and to other high-excitation terms, their departure 
coefficients tend to diverge at log r ^ -2. More importantly, 
similar abundances can be obtained from both the 3905 and 
4102 A lines. 

5. Results 

5.1. Stellar silicon abundances 

The abundance determinations for our program stars are made 
using spectral synthesis. The synthetic spectra are convolved 
with macroturbulence, rotational and instrumental broaden- 
ing profiles, in order to match the observed spectral lines. 
Considering the NLTE abundances, our results do not show a 
large abundance discrepancy between different lines; the final 
abundance scatter of single lines is between 0.01 and 0.14. The 
derived abundances are presented in Table [2] 

5.2. Comparison with other work 

Silicon abundances for metal-poor stars have been determined 
by several groups based on LTE analyses. In Fig|7] we compare 
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Fig. 7. Comparison of derived [Si/Fe] (NLTE) for stars in com- 
mon with other studies. Triangle (A) is from Decauwer et al. 
(2005 ); open circles (o) are from Fulbright d20Q0b ; Filled cir- 
cles (•) are from Gratton et al. (2003); Pluses (+) are from 
Honda et al. ©04); Diamonds (❖) are Ryan et al. (119961 ) 

the [Si/Fe] values (NLTE) determined in this paper with those 
from the literature. Some systematic differences can be seen 
from this figure. In the remaining part of this section we will 
briefly discuss some possible reasons for these differences. 

Decauwer et al. (2005) 

Decauwer et al. (2005) performed a LTE line formation 
analysis for a sample of moderately metal-poor stars, and they 
used weak Si i lines when determining the Si abundances. We 
have one star in common with theirs. Our [Si/Fe] is about 0.1 
dex higher than theirs. The gf values adopted in their study are 
nearly the same as those used here. We note that the optical 
Si i lines are very weak in metal-poor stars, so the determina- 
tion of the continuum is difficult, while the NLTE effects for 
these lines can be negligible. The difference is mostly due to 
the determination of the continuum. 

Fulbright (2000J2002]) 

This analysis deals with a large number of metal-poor 
stars, of which three stars are in common with our sam- 
ple. For the three stars in common with our list we obtain 
A[Si/Fe]=-0.21 + 0.13. We find that the greatest difference 
(~0.35 dex) comes from metal-poor stars such as HD31128. 
They also used weak Si i lines in determining the silicon abun- 
dances, and took the log gf values from the laboratory work 
of Garz (1973). Their gf values are about 0.1 dex lower then 
ours. Also, their effective temperatures are about 200 K lower. 
Thus, the difference in gf values and temperatures can explain 
the discrepancies. 

Gratton etal.(2003) 

The authors determined the silicon abundances for 150 field 
subdwarfs and subgiants. Our results are mostly in agreement 
with theirs. For the six stars in common, the average difference 
is 0.003 + 0.09. 

Honda et al. d2004l > 



Using high resolution, high signal to noise ratio spectra of 
22 very metal-poor stars from Subaru/HDS, the authors con- 
firmed the overabundance ratio [Si/Fe] of metal-poor giants 
found previously by Cayrel et al. (2004). Their results were 
determined from the 4102 A line; the gf value adopted in 
this study is also from the laboratory work of Garz (119731 ). 
It is 0.15 dex lower than ours. We have two stars in common 
with this work; the giant star HD 122563 and the subgiant star 
HD 140283. For the giant, the silicon abundance determined by 
this study is 0.22 dex higher than ours; most of the difference 
comes from the different gf value adopted. For the subgiant, 
their result is -0.04 dex lower than ours; this is due to the large 
NLTE effects (+0. 19 dex) in this star. Thus, the difference in gf 
value is compensated for by the large NLTE effect. 

Ryanetal. (1996) 

This work analysed 19 metal-poor stars using the 3905 A 
line. The gf value adopted in this study was also from the lab- 
oratory work of Garz (119731 ). and it is nearly same as the value 
we adopted. The results of Ryan et al. (1996) are in agreement 
with ours. For the two common stars from Ryan et al. ([1996 ), 
the average difference between theirs and ours is 0.025 + 0.04. 
Part of the difference is due to the NLTE effects for this strong 
line. 



6. Discussion 

6.1. Sii/Sin ionization equilibrium in metal-poor stars 

There are nine stars in our sample with Si abundances deter- 
mined from both Si i and Si n lines, and excellent agreement 
between the two ionization stages is achieved. The average dif- 
ference of A[Sii/Fe]-[Siii/Fe] is 0.001 + 0.04. The largest dif- 
ference between Si i and Si n abundances is found to be 0.09 
for CD -51° 4628. It may due to the large uncertainty in the 
HIPPARCOS parallax for this star (~ 20%). We conclude that, 
within the modelling uncertainties, NLTE leads to consistent Si 
abundances derived from the two ionization stages, while LTE 
fails to give consistent results. 

6.2. Silicon abundance and nucleosynthesis in the 
early Galaxy 

The variation of [Si/Fe] with the stellar metallicity [Fe/H] con- 
tains information about the chemical evolution of the Galaxy. 
Fig. [8] displays the behaviour of the [Si/Fe] ratio (calculated 
in LTE and NLTE) with the metal abundance for all stars con- 
sidered in this paper. One important feature that can be seen 
from Fig. [8] is that there is a large scatter in the LTE results, 
while when the NLTE effects are included, the bulk of [Si/Fe] 
ratio is ~ 0.3 for our program stars. The higher [Si/Fe] ratio 
at lower metallicity found in some studies of metal-poor stars 
(e.g. Fulbright 2002 , Stephens & Boesgaard 2002) is not con- 
firmed by this work. The obvious exceptions are the two moder- 
ately metal-poor stars (CD -51° 4628 and HD 122196), which 
exhibit relatively low [Si/Fe] values. We note that these two 
stars also show lower values of magnesium. Such stars were 
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also found by Nissen & Schuster ( fT997l > and Cohen et al. ( 120071 



Detailed modelling of Galactic chemical evolution has been 
attempted by many authors (e.g. Timmes et al. 119951 Goswami 
& Prantzos 120001 Kobayashi et al. |20061 >. Based on Woosley 
& Weaver's (1995) metallicity-dependent yields, Timmes et 
al. £[995 ) calculate the behaviour of [Si/Fe] as a function of 
metallicity. Their result predicts that [Si/Fe] increases from 

[Fe/H] 0.5 to ~ -1.5, while it decreases from [Fe/H] ~ 

-2 to ~ -3 (see their Fig. 21). The decrease of [Si/Fe] with 
decreasing metallicity is due to (unspecified) uncertainties in 
the low-metallicity massive star models. However, using the 
same yields but with an iron yield reduced by a factor of 2, 
Goswami & Prantzos (2000) show that the [Si/Fe] increases 

from [Fe/H]~ to [Fe/H] 4. They predict that the [Si/Fe] 

ratio is about 0.5 dex at [Fe/H] = -3. Their result suggests an 
increasing [Si/Fe] towards the lower metallicity regime. They 
used a different initial mass function and a different halo model. 
A similar result was found by Kobayashi et al. (2006, see their 
Fig. 10). The Si yields of this work are larger for more massive 
metal-free stars. 



Using magnesium instead of iron as the reference can re- 
move Type la SNe from consideration, because most of the 
silicon and nearly all the magnesium are produced in massive 
stars (Timmes et al. |19951 ). Calculations of nucleosynthesis for 
SNe la show that they should produce some silicon (Nomoto 
et al. 119971 ). The overall behaviour of the [Si/Mg] ratios versus 
[Fe/H] is shown in Fig. |9l where magnesium abundances are 
taken from Gehren et al. (120041 120061 ) and Mashonkina et al. 
(2008 ). Although there is a scatter, our result shows that the 
[Si/Mg] ratio is around -0.1. 



Our observational results provide some implications for the 
nucleosynthesis of silicon. Theoretically, it is expected that sili- 
con is synthesized predominantly in moderate mass (~ 20 Mq) 
Type II SNe, while magnesium is produced primarily in the 
higher mass (~ 35 Mq) Type II SNe (Woosley & Weaver [l995l ). 
The overabundance of silicon and the ratio of [Si/Mg] is about 
-0. 1 in metal-poor stars, suggests that, similar to magnesium, 
silicon is produced by massive type II SNe. The slightly lower 
[Si/Fe] values for extremely metal-poor stars may indicate that 
silicon is produced more by moderately massive type II SNe. 
This result agrees with nucleosynthesis calculations of massive 
stars by Woosley & Weaver ( 119951 ), but it is at variance with 
the calculations by Nomoto et al. ((2006). It should be noted 
that there is only one extremely metal-poor star in our sample; 
it is necessary to investigate more extreme metal-poor stars to 
confirm this suggestion. From a theoretical point of view, since 
silicon is produced in a deeper region than magnesium, the ra- 
tio of Mg/Si is sensitive to the explosion energy and the outer 
border of the mixed region. The higher explosion energy yields 
a smaller Mg/Si ratio, while the Mg/Si ratio decreases with 
the increasing outer boundary of the mixed region (Umeda & 
NomotoEOOSl). 



0.50 E 



0.40 



0.30 r 



33 0.20 r 



o.io : 



o.oo ; 



-3.50 



-3.00 



-2.50 -2.00 
[Fe/H] 



-1.50 -1.00 



Fig. 8. Abundance ratios [Si/Fe] as a function of [Fe/H]. Filled 
circles (•) represent NLTE analysis, while open circles (o) for 
LTE analysis. 
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Fig. 9. Abundance ratios [Si/Mg] for NLTE analysis as a func- 
tion of [Fe/H] . 



6.3. Comparison with DLA abundances 

The chemistry of high-redshift DLA systems should represent 
at least crudely the state of the ISM gas at an early stage in 
the formation of the halo of the Galaxy. It is expected that the 
abundance ratios of elements determined from DLA would be 
consistent with those deduced from spectra of extremely metal- 
poor stars in our Galacy, although the metallicities of DLA sys- 
tems do not reach values as low as those of individual Galactic 
halo stars. The Si/Fe ratio is of crucial importance in under- 
standing the chemical nature of DLAs (Henry & Prochaska 
120071 ). 

As DLA systems are generally in young evolutionary 
stages, it is expected that supersolar [Si/Fe] values will be ob- 
served. Prochaska & Wolfe (2002) found that [Si/Fe] values 
exhibt a plateau 0.3 dex at [Si/H] < -1.5, and Dessauges- 
Zavadsky et al. (2006) found an average [Si/Fe] 0.43 dex for 
their sample, which indicates significant ff-enrichment in the 
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DLA systems at low metallicity. Based on a large number of 
high precision measurements obtained with echelle spectrome- 
ters on 8- to 10-m-class telescopes, Henry & Prochaska (2007) 
also found that [Si/Fe] for all their objects is +0.44. However, 
the mean [Si/Fe] predicted by their model for their sample is 
+0.31 dex, and the ratio is roughly constant. This is very sim- 
ilar to our observational results, but is slightly lower than the 
value obtained from DLA systems. The most likely explana- 
tion is Fe depletion onto dust (Savage & Sembach 1 19961 1 in 
DLA systems. 

7. Conclusions 

We have determined silicon abundances for 14 metal-poor 
stars, spanning the range -3.5 < [Fe/H] <-1.4. All abundances 
are derived from NLTE statistical equilibrium calculations. 
Based on our results we come to the following conclusions: 

1. The [Si/Fe] ratios are overabundant for metal-poor stars, 
and there is an indication that [Si/Fe] slightly decreases 
with decreasing metallicity. 

2. The NLTE effects are different from line to line. The weak 
lines are insensitive to NLTE effects, while the strong 3905 
and 4102 A lines show large NLTE effects. Large depar- 
tures from LTE appear in warm metal-poor stars. The NLTE 
effects increase with increasing temperature, which can 
partly explain the observed results of an anomalous T e g de- 
pendence of [Si/Fe] (Preston et al. 2006; Lai et al. 2008). 

3. A different situation is found for Si n. For the two Si n lines, 
NLTE leads to enhanced absorption in the line cores and 
negative abundance corrections over the range of stellar pa- 
rameters studied here. 

4. The Lya continuum absorption and the radiative bound- 
free cross-sections are very important for Si line formation, 
especially for warm metal-poor stars. 

5. Taking advantage of our statistical equilibrium approach 
and accurate atomic data for the investigated lines, we ob- 
tain good agreement between the Sii and Sin abundances. 
We show that NLTE largely removes obvious discrepancies 
between Si i and Si n obtained under an LTE assumption. 

6. Our results suggest that, similar to magnesium, silicon is 
produced by massive type II SNe. 

It would be important to perform NLTE abundance de- 
terminations for silicon for some warm extremely metal-poor 
stars. 
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Table 2. Stellar silicon LTE and NLTE (for each star, the first and the second row, respectively) abundances given relative to the 
iron LTE abundances derived from the Fe II lines*. 
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Jog g 


IFe/Hl 




3905 


4102 


5690 


5701 


5772 


6142 


6145 


6155 


6237 


6243 


6244 


6347 


6371 


[Si i/Fe] 


[Si n/Fe] 


BD -04° 3208 


6310 


3.98 


-2.23 


1.50 


0.22 


-0.01 
























0.11±0.12 














0.30 


0.24 
























0.27±0.03 




CD -51° 4628 


6076 


4.19 


-1.40 


1 .50 


0.08 


-0.03 




0.09 








0.13 








0.19 


0.18 


0.07±0.05 


0.19±0.01 












0.10 


0.09 




0.09 








0.13 








0.19 


0.18 


0.10±0.01 


0.19±0.01 


HD 29907 


5573 


4.59 


-1.60 


0.90 


0.31 


0.23 




0.32 


0.28 






0.29 


0.31 


0.28 


0.28 


0.31 




0.29±0.02 


0.31 












0.35 


0.29 




0.32 


0.28 






0.29 


0.31 


0.28 


0.28 


0.30 




0.30±0.02 


0.30 


HD31128 


5990 


4.41 


-1.54 


1 .30 


0.34 


0.22 












0.24 


0.24 


0.24 


0.24 


0.24 




0.25±0.03 


0.24 












0.35 


0.34 












0.24 


0.24 


0.24 


0.24 


0.23 




0.28±0.03 


0.23 


HD 34328 


5955 


4.47 


-1.66 


1 .30 


0.37 


0.26 


0.40 










0.34 


0.37 


0.41 


0.39 


0.38 


0.39 


0.36±0.04 


0.39±0.01 












0.42 


0.40 


0.40 










0.34 


0.37 


0.41 


0.39 


0.36 


0.37 


0.39±0.02 


0.37±0.01 


HD 59392 


6020 


3.90 


-1.67 


1.70 


0.24 


0.08 












0.21 


0.22 


0.24 


0.24 


0.26 


0.24 


0.21±0.04 


0.25±0.01 












0.26 


0.24 












0.21 


0.22 


0.24 


0.24 


0.24 


0.23 


0.24±0.01 


0.24±0.01 


tlU /4UUU 


621)3 


4.03 


-2.05 


1 .70 


0.26 


0. 1 9 












0.28 








0.30 




n o/i _i_n (\a 
U.Z4±U.U4 


0.30 












0.30 


0.31 












0.28 








0.28 




0.30±0.01 


0.28 


HD 84937 


6365 


4.00 


-2.15 


1.60 


0.25 


0.11 
























0.18±0.07 














0.30 


0.31 
























0.31±0.01 




HD 97320 


6040 


4.15 


-1.24 


1 .30 


0.32 


0.15 


0.35 


0.32 


0.31 


0.34 


0.33 


0.27 


0.31 


0.31 


0.31 


0.31 


0.34 


0.30±0.03 


0.33±0.02 












0.37 


0.34 


0.35 


0.32 


0.31 


0.34 


0.33 


0.27 


0.31 


0.31 


0.31 


0.30 


0.34 


0.32±0.02 


0.32±0.02 


HD 102200 


6120 


4.17 


-1.28 


1 .50 


0.30 


0.10 


0.34 


0.25 


0.26 


0.27 


0.28 


0.20 


0.20 


0.31 


0.29 


0.25 


0.30 


0.25±0.05 


0.28±0.03 












0.34 


0.28 


0.34 


0.25 


0.26 


0.27 


0.28 


0.20 


0.20 


0.31 


0.29 


0.23 


0.28 


0.27±0.04 


0.28±0.01 


HD 122196 


5957 


3.84 


-1.78 


1.70 


0.12 


-0.04 












0.09 








0.13 


0.13 


0.06±0.06 


0.13 












0.13 


0.10 












0.09 








0.1 1 


0.11 


0.11±0.02 


0.1 1 


HD 122563 


4600 


1 .50 


-2.53 


1.90 


0.24 


0.19 
























0.22±0.05 














0.24 


0.19 
























0.22±0.05 




HD 140283 


5773 


3.66 


-2.38 


1 .50 


0.27 


0.07 
























0.17±0.10 














0.32 


0.26 
























0.29±0.03 




G 64-12 


6407 


4.20 


-3.12 


2.30 


-0.02 


























-0.02 














0.23 


























0.23 





* Our NLTE calculations for Fe I/1I based on the advanced atomic model (Mashonkina et al. 2009 I support the earlier conclusion of Korn et al. i 2003 "i that the NLTE effects for the Fe 
II lines are negligible. 



